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We present a new algorithm for the reduction of one-loop tensor Feynman integrals with n < 4 
external legs to scalar Feynman integrals /„ with n = 3,4 legs in D dimensions, where D = d + 2l 
with integer / > and generic dimension d = 4 — 2e, thus avoiding the appearance of inverse 
t — Gram determinants 04. As long as 04 ^ 0, the integrals I® 4 with D> d may be further expressed 



On 
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by the usual dimensionally regularized scalar functions 7,3 4 . The integrals I® are known at 
04 = 0, so that we may extend the numerics to small, non-vanishing 04 by applying a dimensional 
recurrence relation. A numerical example is worked out. Together with a recursive reduction of 
6- and 5-point functions, derived earlier, the calculational scheme allows a stabilized reduction 
of n-point functions with « < 6 at arbitrary phase space points. The algorithm is worked out 
• ^ explicitely for tensors of rank R<n. 
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1. Introduction 

The efficient and stable evaluation of tensor Feynman integrals, 

with denominators cj, having indices Vy and chords qj, 

Cj = (k-qjf-mj+ie, (1.2) 

is an important ingredient of precision calculations for collider physics. The normalization C(e) 
plays a role for divergent integrals only and is conventional, C(e) = (/^) 2e r(l — 2e)/[T(l +e)r 2 (l — 
e)]. Here, we use the generic dimension d = 4 — 2s and fi = l. 

For n < 4, the problem was basically solved in the Seventies of the last century [1, 2] by tensor 
reduction with an ansatz of linear equations in terms of scalar integrals, and their evaluation in 
terms of logarithms and dilogarithms. For predictions of massive particle production at the LHC 
or ILC, one needs multi-dimensional phase space integrals over typically hundreds to thousands 
of Feynman diagrams with n < 6 external legs and tensor ranks R < n. Further, the so-called 
£-expansion is needed in higher-order calculations. Over the years, a variety of papers appeared, 
and a comprehensive survey is well beyond the scope of this contribution. One of the approaches is 
purely algebraic and rests on the representation of single tensor integrals by scalar integrals, where 
the latter are defined in higher dimensions D = d + 21 and may have also higher indices Vj + m [3], 
e.g.: 

# = /4^n^=-L«r 2 , (i.3) 

J 171 ' r=\ i=l 

d d k " " 1 



# v = / ^^ v n^= L/ijmj^-^c 2 , (i.4) 

^ X = f^n^flc^ = - t q » q ) q ln ijk ^ + \±^q^\ (1.5) 

J 171 r=l uj-k=\ i=l 

= t S?^«^<^ (1.6) 

Tensors like g^ v qf^ are completely symmetrized, and we define = v,y = 1 + <5,y, ra,^ = VijVijki Vijk 
1 + + 8jk etc., and: 

Dstu--- f d^k i^w- 1 

7 M/*-- = / 7^0/2 1 1 l+5„+5 rj +5 rt +--5„-5„-5 r „— ■ • ( L7) 

lli r=l Cr 



The I„ l lY a b ' 5 e -§- * s obtained from ij,^ 1 ' ^ by shrinking line s and raising the powers of inverse 
propagators a,b(s^ a,b); /^'"' },a = . 



2 



Reduction of tensor integrals 



Tord Riemann 



For one-loop tensor integrals, recurrence relations have been derived allowing to represent 
the scalar tensor coefficients l D , , in terms of scalar integrals in generic dimension d and with 
natural indices, usually Vj = 1 [4, 5]. In the present work we need the relation to reduce index j and 
dimension D simultaneously and another relation to reduce dimension D only: 

= -Q 5 /r 2 +i Q n (k-/f- 2 ), d.8) 

0„ {D + l- t V,)J? = Q It 2 - t (°) n (k-/,f- 2 ) . (1-9) 

These relations hold for arbitrary index sets {v,}. The integrals k~/,f and j + /,f are obtained from 1% 
by replacing Vu — > (vt — 1) and Vj — > (v ; - + 1), respectively. For the definitions of Gram determinant 
() n and signed minors like (fy we refer to [6, 7, 8]. If the Gram determinant vanishes, () n = 0, 
relation (1.9) allows to express ^,f{ a | by simpler integrals, some of them with less external legs: 

«{4l() n =o = 7^l(f) ( k "C})" ( L1 °) 

All this works fine, but for tensors of rank R one gets an intermediate scalar basis with dimensions 
up to D = d + 2R, and finally relations in generic dimension d with coefficients (1/()«) R - If () H 
becomes small, and this happens during phase space integrations, numerical instabilities will appear. 

From other approaches it is well-known that the appearance of powers of inverse Gram deter- 
minants for 5-point functions may be avoided completely, and for 4-point functions one has to apply 
special measures if needed. 

In [7], we have demonstrated for pentagons up to R = 3 that a clever use of properties of signed 
minors allows to cancel all the inverse Gram derminants () 5. In [9], we derived a recursive algorithm 
for the representation of (n,R) tensors by (n,R— 1 ) tensors and (n — l,R — 1) tensors, although with 
appearance of inverse Gram determinants ()„. 

Here, we describe an algorithm which combines both approaches and allows an efficient 
evaluation of the tensor integrals in terms of scalar functions; the latter may be evaluated by 
packages like FF [10], LoopTools/FF[l 1, 10], QCDloop/FF [12, 10]. 

The algorithm has been worked out until (n,R) = (6,6) tensors, but it is evident how to go 
beyond that. 

2. A sample reduction free of I/Q5: the 5-point tensor I^ v)l 

For details we have to refer to [8] and references therein. As an example, we consider a rank 
R = 3 pentagon . The rank R = 3 tensor was treated already in [7], but the calculational method 
applied here may be more easily extended to higher ranks. We apply the recurrence derived in [13]: 

where the auxiliary vectors contain inverse Gram determinants: 




(2.2) 



- Us 
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The vector (R = 1) is free of the inverse Gram determinant Q5. Starting the general recursion [9], as 
next the rank R = 2 tensor is expressed by scalar 4-point functions, free of inverse O5: 



MV 



1 







5*=1 



OA rd+2, S ,f0s\ d+2 , s 



+ 



7 ^v 



I_L v (\ i d+2 > s 

'2ffl 5 itW 5 4 



(2.3) 



After further involved manipulations, one may arrive at an expression where also I? is expressed 
by scalar 4-point functions with higher indices and in higher dimensions: 



Eqoj 



£ q ? q y k E Uk +£ g ^ q t ] E QOk , 

fc=l 



(2.4) 



1 



(0)5 * =1 



OA ld +A, 



sk 



+ \° S ) v -/ + 



'J 4,i j [ ' 



(0)5*=! 



2W/ 4 3 U/ 4 



(2.5) 
(2.6) 



The presentation is evidently free of I/Q5, and it is more compact than that given in our earlier 
paper [7]. 



3. The 4-point scalars and their Gram determinants 

We have now to express efficiently the following scalar functions: 

jd+2 jd+A jd jd+2 jd+4 jd jd+2 jd+4 
M >M > y 4,p y 4i >M,i ' 4,ij> l 4,ij > l 4,ij >••• 



(3.1) 



We can again treat only an example. The recurrence relation (1.8) applies e.g. to the scalar function 



4 tit' a PP ear i n g as tensor coefficient in (1.6): 



Cfc)4„ jd+4 , V (k) 
v 'i'\Ji ' 







r(/+4,« , (fe) 4 rrf+4 1 (it) 4 7rf+4 / o n 

a V! i y 3,,7 + 7r y 4,; + 7\ _7 4,i • <A 2 ) 

t=l,t#J (m U4 U4 



Here, we see the appearance of the inverse sub-Gram determinant Q = 04. For every dimensional 
shift, another inverse power of it will appear. In contrast to the case n = 5, this may not be completely 
prevented, but the following strategy is quite useful: Restrict the appearance of 04 to terms related 
to lf +21 , where they may be made implicit, and hold the scalar integrals It,,I2,I\ free of them. 
A lengthy calculation yields: 



(i)0)(°)/j ,u/j , (oj) (°) + (oJQ) + (oO(i) ^ .wrf+4 

ToJToJToyv" - *)«(" + ^4 



o 2 



(°) (d-l)d * /OA /+4 , _ © 4 

s w 3 " h 



Ot 



T d+4,t 



4 ( 0i )( 0t ) + (°j)( 0t ) 
v-> \0k) \0j) ^ \0k) \0i) .d+2,1 J_ 

+ 2- , \ 2 7 3 + /ON 2- 



r=l 



(0)' 



V ; t=i,t#j 



3V 



(3.3) 
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This expression is free of inverse Gram determinants O4. Although, by repeatedly applying (1.9), 
the expression of l£ +6 ,I% +4 in terms of scalar integrals in generic dimension d will introduce 
unavoidably such terms: 

. .,,„„ n ± (L). .,,„„ ,wl 1 

(3.4) 



jd+2l 



\Q) 4 jd+2(l-l) _ y \oUjd+2{l-\),t 

o 4 4 h () 4 3 



d + 2l-3 



The same recurrence relation allows to express these scalar functions I'I +21 in terms of simpler ones 



at 04 = 0, if one writes the recurrence for I 4 
arbitrary 04: 



4+2(1+1) _ 



; see (1.10). One may, however, rewrite (1.9) at 



y 4 



1 



©4 



o 4 (D-3)ir+iQ 



jD,k 



(3.5) 



and apply that relation for small 04. The will be evaluated (approximately) at 04 = 0. This 
should give a better approximation than (1.10) for the higher-dimensional functions at small 04, 
and may be even further iterated. 

In the next section, we will study a numerical example. 



4. Numerical example: Dm 

In order to investigate the stability of the method near a typical kinematical point of vanishing 
sub-Gram determinant, we consider an example given in [14], namely the tensor integrals related to 
a certain box diagram, which in LoopTools [11, 15] notations is: 



D0i(id,0,0,s 3 ,S4,s 12 ,s 23 ,0,M 2 ,0,0). (4-1) 

The Gram determinant is: 

04 = 4 (3) =Det(2^-), (4.2) 

where K[ are the internal momenta, expressible by the (incoming) external momenta pf. K\ = 

p u K 2 = K l + p 2 ,K 3 =K 2 + P3,K 4 = 0. Then, with p] = s t) {p t + pj) 2 = s t j, we get: 

04 = -2*12 [^23 + S3S4 - S 2 3 fa +S4- S U )] . (4.3) 

This Gram determinant vanishes if: 

(523-^3+^12) 

S4 -» S C rU = S23 ~, s • (4.4) 

(^23-^3) 

For 523 = 2 x 10 4 GeV 2 , j 3 = lx 10 4 GeV 2 , s 12 = -4 x 10 4 GeV 2 , it is s crit = -6 x 10 4 GeV 2 , and 
we look at the dependence of the tensor coefficients on 

x=— (4.5) 

Scrit 
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X 




3mD?° 




3m Df n 


0.0 [0] 


-3.15407250453 E- 


10 


-3.31837792633 E- 


10 


- 


- 


10" 


15 [0] 


-3.15407250450 E- 


10 


-3.31837792634 E- 


10 


- 


- 


10" 


10 [0] 


-3.15407189016 E- 


10 


-3.31837792176 E- 


10 


- 


- 


10" 


8 [lin] 


-3.15407250056 E- 


10 


-3.31837790700 E- 


10 


- 


- 


10" 


7 [Hn] 


-3.15407246361 E- 


10 


-3.31837773302 E- 


10 






10" 


6 [lin] 


-3.15407198445 E- 


10 


-3.31837599320 E- 


10 






10" 


5 [lin] 


-3.15405621398 E- 


10 


-3.31835860039 E- 


10 






5 x 


10~ 5 [lin] 


-3.15374472234 E- 


10 


-3.31828141644 E- 


10 






5 x 


10~ 4 


-3.15445603308 E- 


10 


-3.31697733987 E- 


10 


-3.15444393358 E-10 


-3.31694238268 E-10 


7x 


10~ 4 


-3.15403131303 E- 


10 


-3.31683867119 E- 


10 


-3.15408833562 E-10 


-3.31674350685 E-10 


io- 


-3 


-3.15374819252 E- 


10 


-3.31641477212 E- 


10 


-3.15374217567 E-10 


-3.31639655233 E-10 


io- 


-2 


-3.15008000183 E- 


10 


-3.29915920672 E- 


10 


-3.15007998895 E-10 


-3.29915924109 E-10 


10" 


-1 


3.11226750695 E- 


10 


-3.13582331976 E- 


10 


-3.11226750694 E-10 


-3.13582331977 E-10 



Table 1: Dm calculated with the algebraic approach in comparison with LoopTools numerics. Both 
calculations agree in digits shown in boldface. Numbers labeled with [0] rely on (1.10) with 04 = 0, those 
with [lin] rely on (3.5), the others on (3.4). 

In these variables: 

04 = -2^523512(^23-53 + ^12)- (4.6) 

In LoopTools conventions, the tensor coefficients Djji are defined as follows: 

3 3 

iJJ=l i=l 

and for our conventions, see ( 1 .5). The inverse propagators are cj = [{k — qj) 2 — mj] = [(k+Kj- 1 ) 2 — 
rrvj] . Because we assume in our formulae q\ = 0, and in LoopTools it is = 0, one has to care 
about specific correspondences; it is e.g.: 

Din = «222 <J 2 6 2 , (4.8) 

with «222 = V22V222 = 6, and VijVijklfljl given in (3.3). In the example, the tensor coefficients are 
finite. For tensors with non- vanishing l/e" terms, there may arise modifications due to different 
normalizations. 

We present in table 1 and figure 1 the sample numerics for D\ \ \ , for M = 9 1 . 1 876. Our numbers 
for if 222 aie evaluated with a Mathematica notebook and for comparison we also used LoopTools 
v.2.4, both approaches in normal double precision. 

When x < 10~ 3 , the table shows that both the numbers from LoopTools and from a combined 
use of (3.3) and (3.4) (in figure: red line) loose precision significantly. With (3.3) and (1.10) we 
may determine the exact value of Dm at x = 0, and then interpolate to rising x (red broken line). 
The broken blue line uses (3.3) and (3.5), and it extends the region of validity of the broken red line 
considerably, although one does not reach the main region, where the red line is numerically safe. 
A simple numerical interpolation would close the gap quite satisfactory because Dm is a smooth 
function of its arguments near x = 0. With much more effort, namely going further with additional 
iterations a la (3.5), one may reach the same. But, one should have in mind the scales involved, and 
if an accuracy of, e.g., three or four digits is sufficient, simple interpolation schemes will suffice. 
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-3.153x10 



-3.1535x 10 



-3.154 x 10 



-3.1545x10 



-10 



-10 



-10 



-10 - 



-3.155 xlO" 10 L 



-14 -12 -10 



-8 -6 
ln(x) 




Figure 1: The tensor coefficient Dm, defined in (4.7), evaluated with (4.8), as a function of x, defined in 
(4.5). Red line: evaluation of Dm with (3.3), using reductions (3.4); red broken line: use of (1.10); blue 
broken line: use of (3.5). 



5. Summary 

We gave an introduction to our purely algebraic approach to tensor reduction of one-loop 
Feynman integrals. The treatment of inverse sub-Gram determinants has been refined, and a case 
study for vanishing 04 has been presented. For not too small Gram determinants 04, the algorithm 
has been realized in the Fortran package OLOTIC [16], which follows the recursive approach [13] 
for tensor integrals n < 6,R < n. The OLOTIC is being made an open source package. A C++ 
package is under development [17] and will evaluate the tensor integrals according to the scheme 
described here. 
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